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Abstract 

Dynamical quasiparticle properties are determined from lattice QCD along the line 
^ i of the Peshier model for the running strong coupling constant in case of three light 

pH I flavors. By separating time-like and space-like quantities in the number density and 

energy density the effective degrees of freedom in the gluon and quark sector may be 
. specified from the time-like densities. The space-like parts of the energy densities are 

I identified with interaction energy (or potential energy) densities. By using the time- 

■ like parton densities (or scalar densities) as independent degrees of freedom - instead 
. of the temperature T and chemical potential fj,q as Lagrange parameters - variations 
I of the potential energy densities with respect to the time- like gluon and/or fermion 
' densities lead to effective mean-fields for time-like gluons and quarks as well as to 

effective gluon-gluon, quark-gluon and quark-quark (quark-antiquark) interactions, 
i The latter dynamical quantities are found to be approximately independent on 

the quark chemical potential fj,q and thus well suited for an inplementation in off- 
shell parton transport approaches. Results from the dynamical quasiparticle model 
(DQPM) in case of two dynamical light quark fiavors are compared to lattice QCD 
' calculations for the net quark density pq{T,fiq) as well as for the 'back-to-back' 

differential dilepton production rate hy q — q annihilation. The DQPM is found to 
pass the independent tests. 
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1 Introduction 



The 'Big Bang' scenario implies that in the first micro- seconds of the universe the 
entire system has emerged from a partonic system of quarks, antiquaries and gluons 
- a quark-gluon plasma (QGP) - to color neutral hadronic matter consisting of 
interacting hadronic states (and resonances) in which the partonic degrees of freedom 
are confined. The nature of confinement and the dynamics of this phase transition 
has motivated a large community for several decades (cf. [1-3] and Refs. therein). 
Early concepts of the QGP were guided by the idea of a weakly interacting system 
of partons since the entropy density s and energy density e were found in lattice 
QCD to be close to the Stefan Boltzmann (SB) limit for a rclativistic noninteracting 
system [4] . However, experimental observations at the Rclativistic Heavy Ion Collider 
(RHIC) indicated that the new medium created in ultrarelativistic Au+Au collisions 
was interacting more strongly than hadronic matter and consequently this notion had 
to be given up. Moreover, in fine with earlier theoretical studies in Refs. [5-7] the 
medium showed phenomena of an almost perfect liquid of partons [8,9] as extracted 
from the strong radial expansion and elliptic fiow of hadrons as well the scaling of 
the elliptic flow with parton number etc. All the latter collective observables have 
been severely underestimated in conventional string/hadron transport models [10- 
12] whereas hydrodynamical approaches did quite well in describing (at midrapidity) 
the collective properties of the medium generated during the early times for low 
and moderate transverse momenta [13,14]. The question about the constituents and 
properties of this QGP liquid is discussed controversely in the literature (cf. Refs. [15- 
17]) and practically no dynamical concepts are available to describe the dynamical 
freezeout of partons to color neutral hadrons that are finally observed experimentally. 
Since the partonic system appears to interact more strongly than even hadronic 
systems the notation strong QGP (sQGP) has been introduced in order to distinguish 
from the dynamics known from perturbative QCD (pQCD). 

Lattice QCD (IQCD) calculations provide some guidance to the thermodynamic 
properties of the partonic medium close to the transition at a critical temperature 
Tc up to a few times T^., but IQCD calculations for transport coefficients presently 
are not accurate enough [18] to allow for firm conclusions. Furthermore, it is not 
clear whether the partonic system really reaches thermal and chemical equilibrium 
in ultrarelativistic nucleus-nucleus collisions [19] such that nonequilibrium models 
are needed to trace the entire coUision history. The available string/hadron trans- 
port models [20-22] partly fail - as pointed out above - nor do partonic cascade 
simulations [23-26] (propagating massless partons) sufficiently describe the reaction 
dynamics when employing cross sections from perturbative QCD. Some models, e.g. 
the Multiphase Transport Model AMPT [27], employ strong enhancement factors for 
the cross sections, however, use only on-shell massless partons in the partonic phase 
as in Ref . [24] . The same problem comes about in the parton cascade model of Ref . 
[28] where additional 2<-^ 3 processes like gg <-> ggg are incorporated but massless 
partons are considered. 
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On the other hand it is well known that strongly interacting quantum systems re- 
quire descriptions in terms of propagators D with sizeable selfenergies IT for the 
relevant degrees of freedom. Whereas the real part of the selfenergies can be related 
to mean-field potentials, the imaginary parts of 11 provide information about the 
hfetime and/or reaction rate of time-like 'particles' [6]. In principle, off-shell trans- 
port equations are available in the literature [29-31], but have been applied only 
to dynamical problems where the width of the quasiparticles stays moderate with 
respect to the pole mass [32]. On the other hand, the studies of Peshier [33,34] in- 
dicate that the effective degrees of freedom in a partonic phase should have a width 
7 in the order of the pole mass M already slightly above Tg. This opens up the 
problem how to interpret/deal with the space- like part of the distribution functions 
and how to 'propagate' effective degrees in space-time in equilibrium as well as out 
of equilibrium. 

Some elementary steps in answering these questions have been presented in a pre- 
ceeding work [35] where the pure Yang-Mills sector of QCD has been addressed in 
the Dynamical QuasiParticle model (DQPM) of Peshier [33,34]. In the latter work it 
could be demonstrated that the DQPM allows for a rather simple and transparent 
interpretation of QCD thermodynamics and leads to effective strongly interacting 
gluonic degrees of freedom that may be propagated in off-shell transport approaches 
as e.g. in Ref. [36]. The present study is an extension of the work in Ref. [35] to 
three light quark flavors (g = d, s) and to flnite quark chemical potential /ig in 
order to determine the effective partonic degrees of freedom, their mean-fields and 
interactions for a practical implementation in the off-shell parton transport model 
PHSD ^ [36]. First results of the PHSD approach for dynamical phase trajectories in 
central Au-|-Au collisions have been presented in Ref. [37]. 

The outline of the paper is as follows: After a short recapitulation of the dynamical 
quasiparticle model with dynamical quarks in Section 2 novel results on the space- 
like and time-like parts of observables are presented that allow for a transparent 
physical interpretation. In Section 3 we will examine derivatives of the space-like 
part of the quasiparticle energy density with respect to the time-like densities which 
provides information on partonic mean-fields and their effective interaction strength. 
An extension of the model to finite quark chemical potentials is given in Section 4 
including a comparison of the net quark density pg(T,fiq) with IQCD calculations 
in the two flavor case. In Section 5 the 'back-to-back' dilepton production rate is 
calculated in the DQPM and compared to respective IQCD calculations again for 
two dynamical quark fiavors. A summary and short discussion closes this work in 
Section 6. 



Parton-Hadron-String Dynamics 
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2 OfF-shell elements in the DQPM 



2.1 Reminder of the DQPM 



The dynamical quasiparticle model adopted here goes back to Peshier [33,34] and 
starts with the entropy density in the quasiparticle limit [34,38-40], 



-dq 



(lmln(-5-i) + ImE, Re>5,), 



'duj d^p dnp^iyj — jj,q) /T) 
2^(27r)3 df 

■duj d^p dnpiiuj + iia) IT) , , ^ ^ _ „ ^\ 
^ (270^ QT ' (lmln(-5ri) +IniE,Re5,), 



where nsioj /T) = (exp(a;/T) — and nF{{uj — //g)/T) = (exp((tc; — fj,q)/T) + 
denote the Bose and Fermi distribution functions, respectively, while A = (P^ — 11)"^, 
Sq = (P^ — and Sq = (P^ — S^)^^ stand for scalar quasiparticle propagators 

of gluons g, quarks q and antiquarks q. The degeneracy factor for gluons is dg = 
2{N^ — 1) = 16 while for quarks q and antiquarks q we get dq — dq — 2NcNf7/8 
— 15.25 for three flavors Nf. In (1) 11 and T, — T,q T,^ denote the quasiparticle 
selfenergics. In principle, 11 as well as A are Lorentz tensors and should be evaluated 
in a nonperturbative framework. The DQPM treats these degrees of freedom as 
independent scalar fields with a scalar selfenergy 11. In case of the fermions Sq, Sq 
and Eg, Eg (for q and q) have Lorentz scalar and vector contributions but only scalar 
terms are kept in (1) for simplicity which are assumed to be identical for quarks and 
antiquarks. Note that one has to treat quarks and antiquarks separately in (1) as 
their abundance differs at finite quark chemical potential Hq. 

Since the nonperturbative evaluation of the propagators and selfenergies is a formidable 

task (and presently not solved) a more practical procedure is to use physically mo- 
tivated Ansdtze with Lorentzian spectral functions for quarks ^ and gluons. 



E\{uj-Ey + 72 {uj + Ef + 72 




and to fit the few parameters to results from IQCD. With the convention E'^{p) — 
^2 + ^2 — 7^^ the parameters Af ^ and 7 are directly related to the real and imaginary 
parts of the corresponding (retarded) self-energy, e.g. 11 = — 2^70; in case of the 
'scalar' gluons. 

2 In the following the abbreviation is used that 'quarks' denote quarks and antiquarks if 
not specified explicitly. 
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Following [34] the quasiparticle mass (squared) for gluons is assumed to be given by 
the thermal mass in the asymptotic high-momentum regime, i.e. 

M\T) = ^ ((AT, + '-Nf) + f E ^) > (3) 



and for quarks (assuming vanishing constituent masses) as, 



with a running coupling (squared), 

9 (T/Tc) = (11^^ _ 2Nf) ln(A2(T/Te - TjT^y ' ^^-^ 



which permits for an enhancement near Tg [41-43]. Here iVg = 3 stands for the 
number of colors while Nf denotes the number of flavors and /ig the quark chemical 
potentials. The parameters A = 2.42 and Tg/T^ = 0.46 are adopted from [34]. As 
demonstrated in Fig. 1 of Ref. [35] this functional form for the strong coupling 
as = g^/ {4tt) is in accordance with the IQCD calculations of Ref. [44] for the long 
range part of the q — q potential. 

The width for gluons and quarks (for Hg — 0) is adopted in the form [45] 



where c = 14.4 (from Ref. [6]) is related to a magnetic cut-off. Note that in case of 
vanishing number of flavors Nf = the expressions for the masses and the coupling 
reduce to those employed for the pure Yang- Mills sector [34,35]. 

The physical processes contributing to the width 7^ are both gg gg, gq <-> gq 
scattering as well as splitting and fusion reactions gg ^ g, gg ^ ggg, ggg gggg or 
g ^ qq etc. On the fcrmion side elastic fermion-fermion scattering pp ^ pp, where 
p stands for a quark q or antiquark q, fermion-gluon scattering pg ^ pg, gluon 
bremsstrahlung pp ^ pp + g or quark-antiquark fusion qq ^ g etc. emerge. Note, 
however, that the explicit form of (6) is derived for hard two-body scatterings only. 
It is worth to point out that the ratio of the masses to their widths ~ ghi{2c/g'^) 
approaches zero only asymptotically for T — 00 such that the width of the quasi- 
particles is comparable to the pole mass slightly above up to all terrestrial energy 
scales. 

Within the DQPM the real and imaginary parts of the propagators A and S now 
are fixed and the entropy density (1) can be evaluated numerically once the free 
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parameters in (5) are determined. In the following we will assume 3 light quark 
flavors Nf = 3. Since the presently available unquenched IQCD calculations (for 
three flavors) for the entropy density are still accompanied with rather large error 
bars the parameters of the DQPM are taken the same as in the pure Yang-Mills 
sector: A = 2.42, Tg/Tc = 0.46 as determined in Ref. [6]. This is legitimate since 
an approximate scaling of thermodynamic quantities from IQCD is observed when 
dividing by the number of degrees of freedom and scaling by the individual critical 
temperature which is a function of the different number of parton species [46]. 
However, these parameters will have to be refitted once more accurate 'lattice data' 
become available. 

The resulting values for the gluon and quark masses - multiplied by Tc/T - are dis- 
played in Fig. 1 (for fig = 0) by the solid lines while the gluon and quark width 
(7g,7q) - multiplied by Tc/T - are displayed in terms of the dashed lines as a func- 
tion of T/Tc. The actual numbers are found to be quite similar as in the pure Yang 
Mills case [6,33,34] and imply very 'broad' quasiparticles already slightly above T^. 
For fj,q = the ratio 7^/73 =4/9 is the same as for the ratio of the squared masses 
m'^/M'^ = 4/9 and reflects the ratio of the Casimir eigenvalues in color space. Con- 
sequently the ratio of the width to the pole mass is smaller for quarks (antiquarks) 
than for gluons in the whole temperature range. 




Fig. 1. The mass M (solid red line) and width 7^ (dashed blue line) for gluons (l.h.s.) and 
the mass m (solid orange line) and width 7q (dashed green line) for quarks (r.h.s.) as a 
function of T/T^ in the DQPM for A = 2.42, T^/T^ = 0.46, and c = 14.4 (for /Xg = 0). All 
quantities have been multiplied by the dimensionless factor Tc/T. 



In order to fix the scale Tc, which is not specified so far, one may directly address 
unquenched IQCD calculations (for 3 light flavors). However, here the situation is 
presently controversal between different groups (cf. Refs. [47,48] and the discussion 
therein). An alternative way is to calculate the pressure P from the thermodynamical 
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relation, 

dP 



iT • 



by integration of the entropy density s over T, where one may tacitly identify the 
'full' entropy density s with the quasiparticle entropy density s*^^^ (1). Since for 
T < Tc the DQPM entropy density drops to zero (with decreasing T) due to the 
high quasiparticle masses and the width 7 vanishes as well (cf. Fig. 1) the integration 
constant may be assumed to be zero in the DQPM which focusses on the quasiparticle 
properties above T^. 




5 6 7 8 910 



T/T, 



Fig. 2. The DQPM results for e{T/Tc){Tc/T)^ (solid violet line), P{T /Tc){Tc/TY (blue 
dash-dotted line) and the interaction measure W {T /Tc){Tc/T)^ (9) (green dotted line). 
Note the logarithmic scale in T jT^- The energy density e (8) practically coincides with the 
quasiparticle energy density from (16) (dashed red line). 



The energy density e then follows from the thermodynamical relation [41,49] 

e = Ts-P (8) 

and thus is also fixed by the entropy s{T) as well as the interaction measure 

W{T) := e(T) - 3P(T) =Ts-4P (9) 
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that vanishes for massless and noninteracting degrees of freedom. 

The actual resuhs for e- (TjTy are displayed in Fig. 2 (solid violet line), P- {Tc/TY 
(blue dash-dotted line) as well as the interaction measure W ■ {Tc/TY (9) (green 
dotted line) and show the typical pattern from IQCD calculations [46] . The scale 
may now be fixed (estimated) by requiring that the critical energy density e(Tc) is 
roughly the same for the pure Yang-Mills case as for the full theory with dynamical 
quarks in line with the approximate scaling of IQCD [46]. Since e(Tc) in the Yang 
Mills sector is about 1 GeV/fm^ (cf. [35]) the scaled energy density e(T/Tc) from 
Fig. 2 can be employed to fix the critical temperature ~ 0.185 GeV. This leads 
to the 'thumb rule' e ^ 2 (T/TJ^ [GeV/fm^] for T > l.2Tc which is roughly fulfilled 
according to Fig. 2. 




A first test of the DQPM in comparison to IQCD calculations is given for the sound 
velocity (squared), 

= f , (10) 



which does not depend on the absolute (uncertain) scale of Tc. This comparison is 
shown in Fig. 3 by the solid line where the IQCD results have been taken from Ref. 
[50] and correspond to Nf = 2 + 1 with A^^ = 6. The DQPM is seen to reproduce 
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the drop in close to Tc within errorbars and to reach the asymptotic vahic 
= 1/3 approximately for T > 2 Tg. This comparison, however, has to be taken 
with some care since the DQPM assumes massless current quarks whereas the IQCD 
calculations employ finite quark masses. 



2.2 Time-like and space-like quantities 

For the further analysis of the DQPM it is useful to introduce the shorthand notations 
(extending [6,35]) 

• ■ ■ ^ ^^/S ^^''^^'"^ ^^^^ ^^^^/^^ ■ ■ ■ ^^^^ 

with P'^ = — denoting the invariant mass squared. The quark and antiquark 
degrees of freedom, i.e. 2dq — 31.5, are consequently by roughly a factor of two 
more abundant than the gluonic degrees of freedom. The ©(±P^) function in (11) 
separates time-like quantities from space-like quantities and can be inserted for any 
observable of interest. Note, however, that not all space- like quantities have a direct 
physical interpretation. 

We note in passing that the entropy density (1) is dominated by the time-like con- 
tributions for quarks and gluons and shows only minor space-like parts (cf. Refs. 
[34,35]). Furthermore, the entropy density from the DQPM is only 10 - 15 % smaller 
than the Stefan Boltzmann entropy density sg^ for T > 2 Tc as in case of IQCD 
[4]. Since these results provide no novel information an explicit representation is 
discarded. 

Further quantities of interest are the 'quasiparticle densities' 

Nf{T) = trj 1, 7V±(T) = trj 1, N^{T) = trj 1, (12) 

that correspond to the time-like (-I-) and space-like (-) parts of the integrated dis- 
tribution functions. Note that only the time-like integrals over space have a particle 
number interpretation. In QED this corresponds e.g. to time-like photons (7*) which 
are virtual in intermediate processes but may also be seen asymptotically by dilep- 
tons (e.g. e+e" pairs) due to the decay 7* e'^e~{ii^ -\- [20]. 
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Scalar densities for quarks and gluons - only defined in the time-like sector - are 
given by 

Ar;(r) = tr^+ , N^{T)^fr^ (^^y iV|(r) = tr+ (13) 

and have the virtue of being Lorentz invariant. 

Before coming to the actual results for the quantities (12) and (13) it is instructive 
to have a look at the integrand in the quark density (12) which reads as (in spherical 
momentum coordinates with angular degrees of freedom integrated out) 

7(a;,p) = 7^ u;p,{u;,p'') ufHu; - piq)/T) . (14) 



Here the integration is to be taken over u and p from to oo. The integrand I{uj,p) 
is shown in Fig. 4 for T = 1.057; (1-h.s.) and T = 37; (r.h.s.) (/^^ = 0) in terms 
of contour lines spanning both one order of magnitude. For the lower temperature 
the quark mass is about 0.55 GeV and the width 7 k, 0.034 GcV such that the 
quasiparticle properties are close to an on-shell particle. In this case the integrand 
I{u!,p) is essentially located in the time-like sector and the integral over the space- 
like sector is almost neghgible. This situation changes for T — 37; where the mass 
is about 0.7 GeV while the width increases to 7 fa 0.25 GeV. As one observes from 
the r.h.s. of Fig. 4 the maximum of the integrand is shifted towards the line u = p 
and higher momentum due to the increase in temperature by about a factor of three; 
furthermore, the distribution reaches far out in the space- like sector due to the Fermi 
factor npiio/T) which favors small u. Thus the relative importance of the time-hke 
(-I-) part to the space-like (-) part is dominantly controlled by the width 7 - relative 
to the pole mass - which determines the fraction of A*"" with negative invariant mass 
squared (P^ < 0) relative to the time-like part N'^ (P^ > 0). 

The actual results for the different 'densities' (multiplied by (Tc/T)^) are displayed 
in Fig. 5 for gluons (l.h.s.) and quarks (r.h.s.) including the antiquarks. The lower 
(magenta) lines represent the scalar densities A^*, the red solid lines the time-like 
densities A^"*", the green lines the quantities A^~ while the thick solid blue lines are 
the sum A^ = A^+ -|- A^~" as a function of T/T^ (assuming Tg = 0.185 GeV). It is seen 
that A^+ is substantially smaller than A^~ in case of gluons in the temperature range 
1.1 < T/Tc < 10. Except for an overall scale factor the result for gluons is practically 
the same as in the pure Yang-Mills case (cf. Fig. 3 in Ref. [35]). The quantity A^ 
follows closely the Stefan Boltzmann limit Nsb for a massless noninteracting system 
of bosons which is given in Fig. 5 (l.h.s.) by the upper dash-dotted line. Though 
N differs by less than 20% from the Stefan Boltzmann (SB) limit for T > 2Tc the 
physical interpretation is essentially different! Whereas in the SB limit all gluons 
move on the light cone without interactions only a small fraction of gluons can 
be attributed to quasiparticles with density A"+ within the DQPM that propagate 
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p [GeV/c] p [GeV/c] 



Fig. 4. The integrand I{uj,p) (14) as a function of lo and p for quarks at temperatures 
T = 1.05 Tc (l.h.s.) and T = 3 Tc (r.h.s.). At T = 1.05 Tc (l.h.s.) the quasiparticle mass 
amounts to m ~ 0.55 GeV and the width to jq ~ 0.034 GeV while at T = 3 Tc (r.h.s.) 
m w 0.7 GeV and jq ~ 0.25 GeV. The contour hues in both figures extend over one order 
of magnitude. Note that for a convergence of the integrals (12) the upper limits for uj and 
p have to be increased by roughly an order of magnitude compared to the area shown in 
the figure. 

within the lightcone. The space-like part corresponds to 'gluons' exchanged in 
t-channel scattering processes and thus cannot be propagated explicitly in off-shell 
transport approaches without violating causality and/or Lorentz invariance. In case 
of quarks (or antiquarks) the results are qualitatively similar but now the time-like 
part comes closer to the space-like part since the ratio of the width to 
the pole mass {'jq/m) is smaller than the corresponding ratio for gluons as stated 
above. Furthermore, the quantity N is closer to the respective SB limit (for massless 
fermions) due to the lower effective mass of the quarks. In this respect the quarks 
and antiquarks are closer to (but still far from) the massless on-shell quasiparticle 
limit. 

The scalar densities (lower magenta lines) follow smoothly the time-like densities 
A^^ (for gluons as well as quarks+antiquarks) as a function of temperature and 
uniquely relate to the corresponding time-like densities or the temperature T in 
thermal equilibrium. 

The separation of A^"*" and so far has no direct dynamical implications except for 
the fact that only the fraction A^"*" can explicitly be propagated in transport models 
as argued above. Following Ref. [35] we, furthermore, consider the energy densities, 

To± .(T) = Tr± ^ , (15) 
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Fig. 5. The various 'densities' (12) for gluons (l.h.s.) and quarks (r.h.s.) including antiquarks 
(at fiq = 0). The lower magenta lines represent the scalar densities A^'^, the red solid lines 
the time-like densities , the green lines the quantities N~ while the thick solid blue 
lines are the sum = + N" as a function of T/T^. The upper dash-dotted lines 
display the Stefan Boltzmann limits Nsb for reference. All densities are multiplied by the 
dimensionless factor {Tc/T)'^ to divide out the leading scaling with temperature. 

that specify time-like and space-like contributions to the quasiparticle energy densi- 
ties (x = g,q,q)- 
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Fig. 6. The time-like energy density Tqq (red solid line), the space-like energy density Tqq 
(dashed green line) and their sum Tqo = TJq + Tqq (upper blue line) as a function of T/Tc 
for gluons (l.h.s.) and quarks (+ antiquarks) (r.h.s.). All densities are multiplied by the 
dimensionless factor (Tc/T)^. 



The result for the quasiparticle energy densities Tqq and Tqq are displayed in Fig. 6 for 
gluons (l.h.s.) and quarks+antiquarks (r.h.s.) as a function of T /T^. All quantities 
have been multiplied by the dimensionless factor (Tc/T)^ in order to divide out 
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the leading temperature dependence. The lower red solid lines show the time-like 
components Tq{j while the dashed green lines display the space-like parts Tqq which 
dominate over the time-like parts except in the vicinity of T^. The general behaviour 
for the scaled energy densities is similar as for the 'densities' A^='= given in Fig. 5 
since the extra factor cu in the integrand does not change significantly the time-like 
and space-like parts. As in Ref. [35] the space-like parts are interpreted as potential 
energy densities while the time-like fractions are the gluon and quark quasiparticle 
contributions which propagate within the lightcone. 

Summing up the time-like and space-like contributions for gluons, quarks and anti- 
quarks we obtain the total energy density T™, 

— g + ToQ g + T^^g + Tqq_^ -|- T(J^^ -|- Tqq^^ , (16) 



which is displayed in Fig. 2 by the dashed red hne (scaled by Tc/TY). As in case of the 
pure Yang-Mills system in Ref. [35] the quantity T°° practically coincides with the 
energy density e (8) obtained from the thermodynamical relations. Small differences 
of less than 5% show up which indicates that the DQPM in its present formulation 
is not fully consistent in the thermodynamical sense. Since these differences are 
small on an absolute scale and significantly smaller than differences between present 
independent IQCD calculations for 3 quark flavors one may consider T^^{T) ^ e{T) 
and separate the kinetic energy densities TJq from the potential energy densities Tqq 
as a function of the temperature T or - in equilibrium - as a function of the scalar 
densities N'^ or time- like densities iV"*", respectively. 

It is instructive to show the potential energies per degree of freedom Vgg/N^ — 
T^^g/N^ and Ygq/N^ = T^^^/N^ as a function of T/Tc. The corresponding quantities 
are displayed in Fig. 7 (l.h.s.) multiplied by Tc/T in terms of the solid red line 
and the dot-dashed blue line. It is seen that the potential energies per degree of 
freedom steeply rise in the vicinity of T^. and then increase approximately linear with 
temperature T. As expected from the larger width of the gluons the latter also show 
a potential energy per degree of freedom which is roughly a factor of two larger than 
the corresponding quantity for quarks (antiquarks). Consequently rapid changes in 
the temperature (or density) - as in the expansion of the fireball in ultrarelativistic 
nucleus-nucleus collisions - are accompanied by a dramatic change in the potential 
energy density and thus to a violent acceleration of the quasiparticles. It is speculated 
here that the large collective flow of practically all hadrons seen at RHIC [8] might 
be attributed to the early strong partonic forces expected from the DQPM. 

Furthermore, the plasma parameter F defined by the ratio of the average potential 
energy per particle to the average kinetic energy per particle, 

r. = ;^, r, = -^, (17) 

kin,g -'- kin,q 
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Fig. 7. l.h.s.: The potential energies per degree of freedom for gluons Vgg/N^ = Tqq g/iV+ 



(dashed blue line) and for quarks (antiquaries) Vqq/N^ = Tqq q/N^ (solid red line) as a 



function oiT/Tc- All energies are multiplied by the dimensionless factor (Tc/T). r.h.s.: The 
plasma parameter T (17) for gluons (dot-dashed green line) and quarks (solid orange line) 
as a function of T/Tc- Note that F ~ 1 — 2 separates a gas phase from a liquid phase in 
case of Lenard-Jones type of interactions. 



is displayed in the r.h.s. of Fig. 7 for gluons (dot-dashed green line) and quarks (solid 
orange line) as a function of T/Tc. Here the kinetic energy densities are evaluated as 
[6] 



kin,g 



Tr (.; - VP2) 



Tkin,q = Tr {uj - VP2). 



(18) 



The presents results clearly indicate that the plasma parameters F^, F^ are larger 
than unity for both quarks and gluons up to 10 (except for the vicinity of Tc) 
such that the system should be in a liquid phase provided that some attractive 
interaction between the constituents persists. Note that the present evaluation of the 
plasma parameter F is entirely carried out within the DQPM and no longer based on 
estimates for the potential energy as in Refs. [5,6]. Whereas in the earlier estimates 
Tg was dropping with temperature and becoming lower than unity for T > 4 [6] 
the present results indicate than the sQGP should persist for a much larger range in 
temperature (or energy density) and thus also show up in nucleus-nucleus collisions 
at Large Hadron Collider (LHC) energies. Consequently, a partonic liquid is expected 
to be seen also at LHC energies and as a consequence the observed scaling of elliptic 
flow of practically all hadrons with the number of constituent quarks (as seen at 
RHIC) should persist also at LHC. 
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3 Selfenergies and effective interactions of time-like quasiparticles 



Since in transport dynamical approaches there are no thermodynamical Lagrange 
parameters like the inverse temperature (5 = T^^ or the quark chemical potential 
fig, which have to be introduced in thermodynamics in order to specify the average 
values of conserved quantities (or currents in the relativistic sense), derivatives of 
physical quantities with respect to the scalar densities A''^ (or time-like densities 
N^) (x = g,q,q) are considered in the following (cf. Refs. [35,51]). 




Pp [fm"^] Pp [fm"^] 

Fig. 8. The parton potential energy density V (19) - divided by the parton density pp - from 
the DQPM (I.h.s., blue line) as a function of the parton density pp (20). The functional 
form of V{pp) is well reproduced by the expression (22) when divided by pp (green line) 
such that the lines cannot be separated by eye. The r.h.s. shows the gluon fraction a (20) 
from the DQPM (dark blue line) as a function of pp together with the fit (23) (green line). 
Again both lines cannot be distinguished by eye. 

The independent potential energy densities 14 := Tf^ ,^ now can be expressed as 
functions of the scalar densities (or N^) instead of the temperature T (and/or 
quark chemical potential fig). For a determination of mean-field potentials for gluons 
and quarks (antiquarks) it is useful to consider the partonic potential energy density 

V := Too,, + ^00,. + ^00,.- = + Vgg + Vgg (19) 



and to separate a pure gluonic interaction density Vgg from a pure fermionic interac- 
tion density Vgg as well as a gluon-fermion interaction density Vgg. Correspondingly, 
a parton density pp and gluon fraction a is defined via 



p, = N: + N: + Nt, „ = (20) 
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In the present DQPM (for Tc = 0.185 GeV) the parton density pp (20) turns out to 
be a simple function of temperature, 




such that the average distance between the partons is given by d{T /Tc) = Pp^^^ ~ 
(Tc/T)^ °^ [fm] Tc/T [fm]. These relations allow to convert temperatures scales to 
geometrical scales in a simple fashion. 

In Fig. 8 the parton potential energy density V (19) - divided by the parton density 
Pp - is shown as a function of pp (l.h.s.). The functional dependence of V on pp can 
be well approximated by the expression (cf. l.h.s. of Fig. 8) 

V{pp) ^ 0.975pJ-292 _ o.71p;2.i [GeV/fm^] , (22) 



where the numbers in front carry a dimension in order to match the units in GcV/fm^. 
The gluon fraction a (20) is shown on the r.h.s. of Fig. 8 and is well approximated 
by 

a{pp) = 0.29 - 0.075p-°-^^ - 0.15 exp(-1.6pp). (23) 



Adding half of the interaction density Vqp to the gluon part and fermion part sep- 
arately, we have T^^ g — Vgg + Q.bVqg and T^^ g + Tqq ^ = Vgg + 0.5^^^ such that 
Vqg — Vgg =1 = g + g " Tqq g . TYiQ relatlvc fraction of this quantity to the 
total potential energy density then is evaluated as 

/ ^ 

<Pp) = - , y- , y- • (24) 

^ -'- 00,g ^ ^ 00,g ^ 00,3 



Using the Ansatz: 

V99 + % = (1 - 0^ (25) 
then gives 

Vgg^0.5il-C-K)V , Vgg^0.5{l-C + K)V , ^ , (26) 

with still unknown fraction ^ for the interaction density Vgg. 

In order to determine mean-field potentials Ug{pp) for gluons or Uq{pp) for quarks 
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(in the rest frame of the system) one has to consider the derivatives (cf. Ref. [35,52]) 



Ugipp) 



qgJ 



Ugipp 



d{Vgg + V, 



ggJ 



dim + N^) 



(27) 



which by virtue of (25) can be computed as 
19(1 -K + ^)V dpp 



Ugipp 



dpp 



Uqipp) 



ldil + n + i)V 
2 dpp 



dpp 



dim + m^ 



(2J 



The fraction ^ of the interaction density - in principle a function of Pp but here taken 
to be a constant - now can be fixed in comparison to the gluon mean-field from Ref. 
[35] where the pure Yang-Mills sector has be investigated in the same way. This leads 
to ^ ~ 0.3 and separates the total potential energy density V into ~ 26% for the 
gluon-gluon interaction part, 30% for the quark-gluon interaction part (including the 
antiquarks) and ~ 44% for the fermionic interaction part. 
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Fig. 9. The mean-field potentials Ugipp) for gluons (solid red line) and Uqipp) for quarks 
(dashed blue line) as a function of the parton density pp (l.h.s.). The r.h.s. displays their 
ratio as a function of pp which is ~ 2.05 within 10% accuracy. 

The corresponding results for Ugipp) and Uqipp) are displayed in the l.h.s. of Fig. 9 
in terms of the solid red line and dashed blue line, respectively, and show distinct 
minima at pp ~ 2.2 fm~^ which corresponds to an average partonic distance of ~ 0.77 
fm. The actual numerical results for the mean-fields can be fitted by the expressions. 



Ugipp) ^ 70 e-P^'""-^^ + 2.65 pl^^ + 0.45 pj'^ [GeV] 



(29) 



Ugipp) ^ 32 e-''-/°-3i + 1.1 p°-2i + 0.3 pj-^i [GeV] 
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where pp is given in fm~^ and the actual numbers in front carry a dimension in order 
to match to the proper units of GeV for the mean-fields. The ratio Ug/Uq ~ 2.05 as 
can be seen in the r.h.s. of Fig. 9 for a very wide range of parton densities Pp within 
10% accuracy. 

Some comments on the Lorentz structure of the mean fields Ug and Uq appear appro- 
priate. Note that by taking the derivatives with respect to the time-like densities one 
implicitly assumes that a 4-vector current is the physical source of the selfenergies 
and that Ug, Uq are the O'th components of vector fields Ug and Uq {v ~ 0, 1, 2, 3). 
The spatial components are assumed to vanish in the rest frame of the system and 
can be evaluated by a proper Lorentz boost to the frame of interest. This implies 
that the dynamical forces (as space-time derivatives of a Lorentz vector) are Lorentz 
tensors as in case of QED (or vector selfenergies as in the nuclear physics context 
[52]). On the other hand one might consider the notion of purely scalar selfenergies 
where derivatives of the potential energy density with respect to the scalar density 
(e.g. dT(^ ,j./dN^) define an effective mass M* [52]. In a selfconsistent framework 
then the quasiparticles masses (4) should be given by the derivatives with respect 
to the scalar densities. However, this relation is not fulfilled at all since a numerical 
evaluation of the scalar derivatives gives effective masses that are larger by more 
than an order of magnitude as the quasiparticle masses introduced in (4)! This re- 
sult might have been anticipated since the effective forces for a gauge (vector) field 
theory should be dominated by Lorentz forces as in case of QED. For completeness 
we give the fitted expressions for the quasiparticle masses (4) as a function of the 
scalar densities, respectively, 

Mg{N'g) ^ 0.7{N'gf-^^^ + 0.102(Ar^^)-2/3 [GeV] , (30) 



mAN^^q) ^ OMiN^^X"'' + 0-14(A^,%-)-^/^ [GeV], 

where the scalar densities are given in units of fm^^ and the numbers in front carry 
a dimension to match the units of GeV for the masses (using -/V^+g — + N?). 

Some information on the properties of the effective gluon-gluon, quark-gluon and 
quark-quark interaction may be extracted from the second derivatives of the potential 
energy density V, i.e. 

V (o)- ^'^^^ - l d\l-i-n)V ( dp, y 

^..(Pp)-^~2 d^l \dNj) ' ^^^^ 

d^q \d\\-i^K)V I dpp V 



VqqiPp) 
Vqg{pp) :-- 



d\m ( dpp \ ( dpp 
a(iV+ + N^)dN+ dpi [d{N+ + N^)) [dN^ 
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The numerical results for the interactions (31) are displayed in Fig. 10 (l.h.s.) for 
the effective gluon-gluon (solid red line), gluon-quark (solid blue line) and quark- 
quark interaction (dashed green line). All interactions show up to become strongly 
attractive at low parton density Pp < 2.2 fm~^, change sign and become repulsive 
for all higher parton densities. This situation has been the same in the pure Yang- 
Mills case [35] (at slightly lower gluon density ^1.4 fm~^). Note that the change 
of quasiparticle momenta (apart from collisions) will be essentially driven by the 
(negative) space-derivatives —VUj{x) = —dUj{pp)/dpp V pp{x) which implies that 
the partonic quasiparticles (at low parton density) will bind with decreasing density, 
i.e. form 'glueballs', mesons, baryons or antibaryons dynamically close to the phase 
boundary and repel each other for pp > 2.2 fm~'^. Note that color neutrality is 
imposed by color-current conservation and only acts as a boundary condition for the 
quantum numbers of the bound/resonant states in color space. 




1 10 100 1000 10 100 1000 

pp [fm"'] pp [fm"'] 

Fig. 10. l.h.s. The effective gluon-gluon (solid red line), gluon-quark (solid blue line) and 
quark-quark interaction (dashed green line) from the DQPM for ^ = 0.3 (see text). The 
r.h.s. displays the ratios Vgg/vqq (solid red line) and Vqg/vqq (solid blue line). The straight 
light blue line is the ratio of the Casimir eigenvalues, i.e. Cg/Cq= 9/4. 

The r.h.s. of Fig. 10 displays the ratios Vgg/vqq (solid red line) and Vqg/vqq (solid blue 
hne) as a function of the parton density pp and demonstrates that Vqg/vqq ~ 9/4, 
which is the ratio of the Casimir eigenvalues. 

A straight forward way to model the parton condensation or clustering to confined 
glueballs or hadrons dynamically (close to the phase transition) is to adopt screened 
Coulomb-like potentials Vc{r, A) with the strength / (Pr Vc{r, A) fixed by the interac- 
tions Vgg{pp) , Vqg{pp) , Vqq{pp) from (31) and a screening length A from IQCD studies. 
For the 'dilute parton regime' {pp < 2.2 fm~^), where two-body interactions should 
dominate, one may solve a Schrodinger (Dirac or Klein-Gordon) equation for the 
bound and/or resonant states. This task is not addressed further in the present 
study since for the actual applications in the PHSD approach [36] the formation of 
glueballs is discarded and the formation of resonant hadronic states close to the phase 
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boundary is described by density- dependent transition matrix elements (~ \vqq{pp)\^) 
between partons of 'opposite' color with fixed flavor content. In this way the energy- 
momentum conservation, the flavor current conservation as well as 'color neutrality' 
are exphcitly fulfilled in the (PHSD) transport calculations for interacting particles 
with spectral functions of finite width. 



4 Finite queirk chemical potential /x. 

The extension of the DQPM to finite quark chemical potential fig is more delicate 
since a guidance by IQCD is presently very limited. In the simple quasiparticle model 
one may use the stationarity of the thermodynamic potential with respect to self- 
energies and (by employing Maxwell relations) derive a partial differential equation 
for the coupling g'^{T,fiq) which may be solved with a suitable boundary condition 
for g'^{T,fiq = 0) [41]. Once g'^{T,fiq) is known one can evaluate the changes in 
the quasiparticle masses with respect to T and fiq, i.e. dM^/d/iq and dM^/dT (for 
X — QiQi Q) and calculate the change in the 'bag pressure' AS (cf. Refs. [41,53] for 
details). However, such a strategy cannot be taken over directly since additionally 
the quasiparticle widths 7^(7^, A*^) have to be known in the (T, //g) plane in case of 
the DQPM. 

In hard-thermal-loop (HTL) approaches [54,55] the damping of a hard quark (or 
gluon) does not depend on the quark chemical potential explicitly [57] and one might 
employ (6) also at finite /x. This, however, has to be considered with care since HTL 
approaches assume small couplings g"^ and should be applied at sufficiently high 
temperature, only. Present IQCD calculations suggest that the ratio of pressure to 
energy density, P/ e, is approximately independent on ^q as a function of the energy 
density e [56]. Accordingly, the functional dependence of the quasiparticle width 7 
on jiq and T has to be modeled in line with 'lattice phenomenology' (see below). 

Assuming three hght flavors {q = u, d, s) and all chemical potentials to be equal 
{/lu = IJ'd = l^'s = lA equations (3) and (4) demonstrate that the effective gluon and 
quark masses are a function of 



Since the coupling (squared) (5) is a function of T/T,, a straight forward extension 
of the DQPM to flnite 11 is to consider the coupling as a function of T* /Tc{ij) with 
a //-dependent critical temperature. 




(32) 



(/. = 0)(1 



_1 ^ 

27r2 r,(0) 




(33) 
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The coefficient in front of the //^-dependent part can be compared to IQCD cal- 
culations at finite (but small) /i which gives 0.07(3) [58] instead of 0.05 in (33). 
Consequently one has to expect an approximate scaling of the DQPM results if the 
partonic width is assumed to have the form (6), 

j(Ta)^N 9\TVU^)) ^ 2c 



,^ , A^c - 1 9^{T*lTAp)) _ , 2c 



2iVc 47r g\T*lT,{ii)) ' 

where g^{T jT^) has been replaced by g'^{T* /Tc{ix)). In fact, as will be demonstrated 
below, this choice leads to an approximate independence of the potential energies 
per degree of freedom as a function of /ig. Nevertheless, the conjecture (34) should 
be explicitly controlled by IQCD studies for Nf=3 at finite quark chemical poten- 
tial. Unfortunately, this task is presently out of reach and one has to live with the 
uncertainty in (34) which is assumed in the following investigations. 

Within the scaling hypothesis (33), (34) the results for the masses and widths in 
Section 2.1 stay about the same as a function of T*/Tc{fi) when dividing by the tem- 
perature T. This also holds approximately when displaying the masses and widths as 
a function of the parton density pp for different chemical potentials /i as demonstrated 
in Fig. 11. The latter quantities can well be fitted by the expressions 

Mgipp) ^ 0.41 + 0.38 [GeV] , (35) 



7,(pp) ~ 0.235 pi''' - 0.14 p;' [GeV] , 

2 4 
rriqiPp) ~ -^Mgipp), -fq{pp) ^ - -fg{pp), 

with Pp given in units of fm~^. Note that according to the parametrization (35) the 
width might become negative for very small Pp] in actual transport applications it 
should be set to zero in such cases [36]. 

The more interesting question is how the energy density e (8) and the pressure 
P (from (7)) change with quark chemical potential p — pq in the DQPM. This 
information is provided in Fig. 12 where the upper l.h.s. shows the energy density e 
(8) (scaled in terms of Tco = Tc{p = 0) = 0.185 GeV ) as a function of T* /Tc{p). Here 
a scaling of the 'temperature' T* with Tc(/i) (33) is used since the phase boundary 
changes with the quark chemical potential p. The energy density e is seen to scale 
well with (T/TcoY as a function of temperature for T*/Tc{p) > 3, however, increases 
slightly with p close to the phase boundary where the scaling is violated on the level 
of 20%. This violation in the scaling (seen in the upper left part of the figure) is 
essentially due to an increase of the pressure P which is displayed in the lower left 
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Fig. 11. The gluon mass M and width 7^ (l.h.s.) as a function of the parton density pp for 
various chemical potentials fiq = n from fi = to n = 0.21 GeV in steps of 0.03 GeV. The 
r.h.s. displays the mass m for quarks and width 7^ for the same quark chemical potentials. 

part of the figure as a function of T*/Tc(/i) for the same chemical potentials /i from 
fi = to 0.21 GeV in steps of 0.03 GeV. Note that a quark chemical potential of 
0.21 GeV corresponds to a baryon chemical potential of fiB = 3/i =0.63 GeV which 
is already substantial and the validity of (33) becomes questionable. 

Since the pressure P is obtained from an integration of the entropy density s over 
temperature (7) the increase in P with /i can directly be traced back to a correspond- 
ing increase in entropy density. The latter is dominated by the time-like quasiparticle 
contributions thus 'counting' the effective degrees of freedom, 

= iV; + iV+ , , iV;+, = + iV+ . (36) 



The upper r.h.s. of Fig. 12 shows pp versus T*/Tc{p) (multiplied by T^q/T'^) for 
chemical potentials from /i = to /i = 0.21 GeV in steps of 0.03 GeV. Indeed the 
scaled parton density increases with /i; this increase is most pronounced for lower 
temperatures and becomes substantial for /i = 0.21 GeV. One should recall, however, 
that a tri-critical endpoint (in the T, /x plane) is expected for /i^ = 3/i ~ 0.4 GeV 
[59] which corresponds to p ~ 0.13 GeV. When restricting to the interval < /x < 
0.13 GeV the explicit change in the parton density Pp with p stays very moderate. 
This also holds for the energy density e and the pressure P. 

Quite remarkably the potential energy per time-like fermion T^Q ^^^/N^j^q changes 
very little with p as can be seen from the lower right part of Fig. 12 where the 
latter quantity is displayed for the same chemical potentials as before (lower green 
lines) as a function of T*/Tc{p). This also holds for the potential energy per time-like 
gluon T(^ g/N^ (upper magenta lines). Accordingly the potential energy per time-like 
degree of freedom is essentially a function of T*/Tc{p) alone. 
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Fig. 12. The energy density e{T* /Tc{^)) (8) for quark chemical potentials /i from to 
0.21 GeV in steps of 0.03 GeV as a function of the scaled temperature T* /Tc{ij,) for Nf 
=3 (upper l.h.s.). The pressure P(T* /Tc{fi)) from the thermodynamical relation (7) for 
quark chemical potentials from to 0.21 GeV in steps of 0.03 GeV as a function of the 
scaled temperature T* /Tc{fJ,) for Nf =3 (lower l.h.s.). The parton density pp (for quark 
chemical potentials fi from to 0.21 GeV) is shown in the upper r.h.s. The potential 
energy per time-like 'gluon' (upper magenta lines) and the potential energy per time-like 
'quark+antiquark' (lower green lines) for the same quark chemical potentials are displayed 
in the lower r.h.s. Note that e and P are scaled by the dimensionless factor {Tco/T)^ (with 
Tco = 0.185 GeV) while the parton density pp is scaled by the factor (Tco/T)^. 

The time-like densities for the fermions N^^^^ (36) (upper l.h.s., dark blue lines) and 
the space-like quantities N~j^^ (upper r.h.s., light blue lines) are shown in Fig. 13 as a 
function of T*/Tc{fi) for chemical potentials fi from /i = to /i = 0.21 GeV in steps 
of 0.03 GeV (scaled by T^q/T^). Both quantities increase with yU in a comparable 
fashion such that their ratio stays approximately constant for T* > 1.5 Tc{fi). The 
space-like quantity A*"" (lower l.h.s., dashed green lines) practically is independent 
from fi as well as the time-like gluon density A^"^ (lower l.h.s., red lines). Accordingly 
the gluon fraction 



a(T,/i) 



N+ 

9 



(37) 



q+q 



which is displayed on the lower r.h.s. of Fig. 13 for the same set of quark chemical 
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Fig. 13. Upper parts: The time-like fermion density -/Vg+g (l.h.s., dark blue lines) and the 
space-like quantity N~j^g (r.h.s., light blue lines) for quark chemical potentials fi from to 
0.21 GeV in steps of 0.03 GeV as a function of the scaled temperature T* /Tc{n). Lower 
parts: the time-like gluon density (l.h.s., lower red lines) and the space-like quantity 
N~ (l.h.s., dashed green lines) for the same quark chemical potentials as a function of 
T*/Tc{iJ.). Both, N~ and change only very little with in the whole temperature 
range; the relative changes are of the order of the line width. The gluon fraction (37) is 
shown on the r.h.s. as a function of the scaled temperature T*/Tc(/x) for quark chemical 
potentials from to 0.21 GeV in steps of 0.03 GeV. Note that all 'densities' have been 
multiplied by the dimensionless factor (Tco/T)'^. 

potentials as a function of T*/Tc{fi), decreases with fi. It drops to zero below the 
phase boundary because the difference between the gluon effective mass and the 
fermion effective mass becomes large below Tc. 



We continue with the time-like and space-like energy densities for the fermions and 
gluons as a function of fi and T*/Tc{fi) which are displayed on the l.h.s. of Fig. 14. The 
upper part shows that the space-like and time-like energy densities for the fermions 
increase with fi roughly in a similar fashion as the fermion 'densities' such that their 
ratio is approximately independent on /i. The lower part of Fig. 14 (l.h.s.) indicates 
that the space-like energy density for gluons (green lines) as well as the time-like 
energy density for gluons (red lines) are approximately independent from /i within 
line width. When separating the time-like fermion energy density into contributions 
from quarks q and antiquarks q (r.h.s. of Fig. 14) we find an increase of ^ with /i 
which is not fully compensated by a decrease of Tgjj ^ with quark chemical potential. 
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Fig. 14. l.h.s.: The time-like fermion energy density (blue lines) and the space-like 

fermion energy density Tq^ ^_|_^ (dashed red hnes) for quark chemical potentials /x from 
to 0.21 GeV in steps of 0.03 GeV as a function of the scaled temperature T* /Tc^n) (upper 
l.h.s.). The lower l.h.s. shows the time-like (dashed red lines) and space-like gluon energy 
density (green lines) for the same chemical potentials as a function of the scaled tempera- 
ture. Both, Tqq g and T^^ ^ are roughly independent from n in the whole temperature range. 
The time-like (upper r.h.s) and space-like energy densities (lower r.h.s.) for quarks (upper 
red lines) and antiquarks (lower green lines) are separately displayed for quark chemical 
potentials /i from to 0.21 GeV in steps of 0.03 GeV for Nf = 3. All quantities have been 
multiplied by the dimensionless factor T^^q/T^. 

Similar dependences on /i and T* are found for the space-like sector (lower part, 
r.h.s.). Thus when summing up all space-like and time-like energy densities from 
the fermions and the gluons a small net increase in the total energy density with fi 
survives (see also upper left part of Fig. 12). 



We note in passing that derivatives of the various energy densities w.r.t. the time-like 
gluon or fermion densities (as investigated in detail in Section 3) are approximately 
independent of fi such that the effective potentials Ug{pp), Ug{pp) and Ug{pp) (28) 
stay practically the same. Since this result may be inferred already from the /i-(in)- 
dependence of the potential energy per time-like degree of freedom (displayed in the 
lower r.h.s. of Fig. 12) an explicit representation is discarded. This implies that the 
mean-fields (28) or parametrizations (29) may be employed also at finite (moderate) 
net quark density A^^'" — A*"^ which simplifies an implementation in parton transport 
models (as e.g. PHSD). 
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Whereas for vanishing quark chemical potential /i=0 the quark and antiquark den- 
sities are the same this no longer holds for finite fi where the differences 

pf = - (38) 



are of separate interest since pq = Pg + Pq is the zero'th component of a conserved 
flavor current (separately for each flavor u,d,s, ..).ln order to obtain some idea about 
space-like and time-like contributions of the net quark density pq in the DQPM we 
first plot the time-like component in the l.h.s. of Fig. 15 as a function of T*/Tc{p) 
for different p from to 0.3 GeV (as before). In order to divide out the leading 
scaling with temperature the time-like densities have been multiplied by Tj^/T^ on 
the l.h.s. of Fig. 15 (as known from Fermi systems in the nuclear physics context). 
The actual results show an approximately linear increase in p (l.h.s.) which suggests 
to study the scaled quantities p^ ■ T^q/^pT"^). However, a better scaling is obtained 
when multiplying p^ by T'^^/ [pT*"^). The latter quantities are displayed in the r.h.s. 
of Fig. 15 for the same set of chemical potentials as before. This approximate scaling 
allows to estimate the net quark density as 

^ 10 /i T*VTi, [fm-3] (39) 




Fig. 15. l.h.s.: The time- like quark density Pq (38) (scaled by the dimensionless factor 
T^q/T"^) for quark chemical potentials fi from to 0.21 GeV in steps of 0.03 GeV as a 
function of the scaled temperature T*/Tc{p.). r.h.s.: The time-like quark density (38) 
(red solid lines) and the space- like quantity p~ (38) (dashed green lines) for the same quark 
chemical potentials p. In this part of the figure all quantities are scaled by the factor 

The approximate scaling depicted in Fig. 15 is a prediction of the DQPM in case of 
3 quark flavors and should be controlled by IQCD calculations. For two light flavors 
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{Nf = 2) some comparison can be made in the low temperature (and low fi) range. 
Unfortunately the simple scaling relation (39) does not hold at small T such that 
an explicit comparison has to be presented between the DQPM for Nj = 2 and the 
IQCD calculations from Ref. [60]. The latter calculations have been carried out on 
a 16^ X 4 lattice with two continuum flavors (of p4-improved staggered fermions) 
with mass m = OAT. Though the fermion masses are not really 'light' in the IQCD 
calculations the actual IQCD results may serve a test of the present DQPM in the 
2- flavor sector (for low fi and T). The explicit IQCD results for the net quark density 
Pq (divided by T^) are displayed in Fig. 16 as a function of T /T^o in terms of the 
various symbols. The symbols of equal color correspond to p/Tco = 1.0, 0.8, 0.6, 0.4, 
and 0.2 from top to bottom and by eye show an approximately linear dependence 
on fi. The explicit results from the DQPM for Nf = 2 are presented in terms of the 
dark green lines and approximately follow the IQCD results (at least for smaller fi). 
Since the systematic errors of the IQCD calculations are not known to the author 
an explicit refitting of the parameters X^T^/T^ and c for Nf =2 (cf. Section 2.1) 
is discarded here. Nevertheless, the qualitative (and partly quantitative) agreement 
between the DQPM and IQCD provides a test for the basic concepts of the DQPM. 




Fig. 16. The quark density pq (38) (divided by the factor T^) for quark chemical potentials 
/x/Tco from to 1.0 in steps of 0.2 as a function of the scaled temperature T/Tcq from the 
DQPM (green lines) for Nf =2. The different symbols represent the IQCD results from Ref. 
[60] for the same quark chemical potentials as a function of T/Tcq (also for Nf =2). Note 
that the temperature axis here is given by T/Tco = T/Tc{fj, = 0) and not by T* /Tc{p)\ 
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5 Dilepton radiation from the sQGP 



The properties of the sQGP so far have been fixed in the DQPM by specifying 
the (vector) selfenergies/potentials as well as effective interactions for the time-like 
partons. As shown in Ref. [6] this leads to a strongly interacting partonic system with 
a shear viscosity to entropy density ratio close to r^/s ~ 0.2. However, the predictions 
from the DQPM should be controlled by independent IQCD studies to get some idea 
about the reliability of the approach. As mentioned before transport coefficients like 
the shear viscosity t] are available from IQCD [18] but the present accuracy is not 
satisfactory. On the other hand some IQCD information is available from the Bielefeld 
group on the electromagnetic correlator which is initimately related to the dilepton 
emission rate [61-63]. 
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Fig. 17. The 'back-to-back' dilepton emission rate (41) from the DQPM (solid red line) in 
comparison to the Born approximation (40) for massless partons (dashed blue line) and 
the results from the IQCD analysis in Ref. [63] at T = 1.5 Tc in case of Nf = 2 (full dots). 

In order to provide some further comparison to IQCD results we consider the dilepton 
production rate in thermal equilibrium (at temperature T) in a two-ffavor QGP as 
in Ref. [63]. For massless quarks and antiquarks the emission rate of 'back-to-back' 
leptons (e+e~) is given by [63,64] 



dW 



ha u u , 

= °) = 36^ ^^^2t) "^^2T' 



(40) 



where oj is the invariant mass of the lepton pair and np denotes the Fermi distribution 
function. In (40) a is the electromagnetic coupling constant. For massless partons, 
furthermore, the magnitude of the lepton momenta is given by |p| = uj/2 while their 
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direction is opposite in the dilepton rest frame. Neglecting the rest mass of leptons 
their energy is u;/2 in the dilepton rest frame, too. The expression (40) changes in 
case of spectral functions with finite width to 



dW 5a^ 



oo CXJ CXJ 

j dwi j duj2 j dp —— pq{uJi,p;T)pq{uj2,p;T) (41) 





5{u; -ui- u;2) nF{—) ^^(777) 



^2 V ryrp; ryrp. 

with P? —cUj—p'^ denoting the invariant mass of the annihilating partons j — 1,2. 

In (41) the factor \/\/uj'^ gives a flux correction in case of massive quasiparticles with 
A(a;, y,z) ^ {x - y - zf - Ayz. 

Since the spectral functions in the DQPM are fixed (cf. Section 2) the lepton emission 
rate (41) can be evaluated without introducing any further assumption (or parame- 
ter). The results for the different emission rate for Nf = 2 arc shown in Fig. 17 for 
T = 1.5 Tc (solid red line) in comparison to the limit (40) (dashed blue line) and the 
IQCD results from [63] (full dots). The IQCD dilepton rate has been obtained from 
the temporal correlators in IQCD employing the 'maximum entropy method' which 
has a systematic error in the order of 30 to 50 % [63] depending on the energy scale 
considered. The results from Fig. 17 demonstrate a drastic suppression of low mass 
lepton pairs due to the finite mass of the partons. On the other hand the spectra 
from the DQPM are in qualitative agreement with the IQCD results from Ref. [63] 
when including the systematic uncertainties of the latter approach. Thus the DQPM 
(in its present version) passes a further test in comparison to IQCD. 



6 Conclusions and discussion 



The present study has provided a novel interpretation of the dynamical quasipar- 
ticle model (DQPM) by separating time-like and space-like quantities for particle 
densities, energy densities, entropy densities ect. that also paves the way for an off- 
shell transport approach [36]. As known from previous studies [34,35] the entropy 
density s is found to be dominated by the on-shell quasiparticle contribution while 
the space-like part of the off-shell contribution gives only a small (but important) 
enhancement. However, in case of the parton 'densities' A^^; = + and the 
energy densities Too,^ = Tqq -|- Tqq^ (,t = g, q, q) the situation is opposite: here the 
space-like parts {N^,Tqq^^) dominate over the time-like parts {N~^,T(^) except close 
to Tc where the independent quasiparticle limit of on-shell particles is approximately 
regained. The latter limit is a direct consequence of the infrared enhancement of the 
coupling (5) close to Tc (in line with the IQCD studies in Ref. [44] ) and a decrease 
of the width 7 (6) when approaching Tc from above. 
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Since only the time-like part of the parton density can be propagated within the 
lightcone the space-like part A^^ has to be attributed to t-channel exchange partons 
in scattering processes that contribute also to the space-like energy densities T^^,,.. 
The latter quantities may be regarded as potential energy densities V^. This, in fact, 
is legitimate since the total quasiparticle energy density (16) very well matches 
the energy density (8) obtained from the thermodynamical relations. Only small 
deviations indicate that the DQPM in its straightforward application is not fully 
consistent in the thermodynamical sense. 

By taking derivatives of the potential energy densities T4 with respect to the time- 
like gluon and fermion densities one may deduce mean-field potentials Ux for gluons 
and quarks (antiquarks) as a function of the parton density pp which enters instead 
of the thermodynamical Lagrange parameters T and fig. Second derivatives w.r.t. 
the gluon and/or fermion densities then define effective interactions between gluons, 
quarks and quarks and gluons. We find that the gluon-gluon interaction is stronger 
than the quark-quark (antiquark) by roughly a factor of 4-5 whereas the quark-gluon 
interaction strength approximately scales (relative to the quark-quark interaction) 
by the ratio of Casimir eigenvalues, i.e. 9/4. The effective parton-parton interactions 
are repulsive for pp > 2.2 fm~^ and turn strongly attractive at lower density. Since 
at low parton densities (or average distances larger than 0.77 fm) the system may 
be described by 2- and 3-body dynamics the strong attraction will lead to color 
neutral bound states of gluons (glueballs), quark and antiquarks (mesons) as well as 
3 quarks or 3 antiquarks (baryons or antibaryons). The mean-field potentials Ux are 
found to be independent on the quark chemical potentials Pg (within 5%) such that 
the respective analytical approximations (29) are well suited for an implementation 
in off-shell parton transport approaches (as e.g. PHSD). 

Since for the case of three light flavors (iV/=3) the effective masses (squared) in 
the DQPM scale with T*^ := + /U^/vr^ and the spectral width of the dynamical 
quasiparticles is assumed to be given by the 'modified' HTL expression (34) an 
approximate scaling of various quantities is found in the DQPM at finite quark 
chemical potential Pg. As an example we find that the potential energy per time- 
like gluon as well as the potential energy per time-like fermion are approximately 
independent on pg which directly leads to the independence of the mean-fields Ux 
on Pg (as mentioned above). Furthermore the net quark density is found to roughly 
scale as Pg ~ PgT*'^ which might be directly checked by IQCD. In addition the energy 
density e and the pressure P are found to increase only slightly with pg close to T(. 
such that the ratio P/e (as a function of e) is also approximately independent on pg. 

Explicit comparisons of the DQPM calculations with IQCD results have been per- 
formed for the sound velocity (squared) (10) in the 2+1 fiavor case and in the two 
flavor sector {Nf —2) for the net quark density as a function of temperature T and Pg 
from [60]. The results are in an acceptable agreement such that the DQPM - which 
presently has been flxed only in the pure Yang-Mills sector {Nf = 0) - reproduces the 
dominant dependences on T and pg. Furthermore, it could be shown that the 'back- 
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to-back' dilepton emission rate from the DQPM in case of Nf—2 is in qualitative 
agreement with IQCD studies from Ref. [63] at T = 1.5 Tc thus providing a further 
independent test of the modeL Note that in a conventional quasiparticle model with 
vanishing width jg the dilepton emission rate only gives contributions for invariant 
masses above 2 ruq which is not in line with the IQCD result from Ref. [63]. 

Some note of caution with respect to the present DQPM appears appropriate: the 
parameters in the effective coupling (5) and the width (6) have been fixed in the 
DQPM by the entropy (1) to IQCD results for Nf—0 assuming the form (2) for the 
spectral function p{uj). Alternative assumptions for will lead to slightly different 
results for the time-like and space-like densities, energy densities etc. but not to a 
qualitatively different picture. Also it is presently unclear if the three parameters 
(A, Ts/Tc, c) employed in Section 2.1 are approximately the same in case of two or 
three dynamical flavors. More precise calculations from IQCD should allow to put 
further constraints on the form of the spectral function p{uj) and to fix the basic 
model parameters in the effective coupling. Also the 'modified' HTL expression (34) 
for the quasiparticle width has to be controlled by IQCD at finite pg as well as 
transport coefficients like the shear viscosity rj or related correlators. 

The author acknowledges valuable discussions with E. L. Bratkovskaya, M. H. Thoma 
and A. Peshier. 
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